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Abstract 

We consider several statistical models defined on the Farey fractions. Two of these models may be 
regarded as "spin chains" , with long-range interactions, while another arises in the study of multifractals 
associated with chaotic maps exhibiting intermittency. We prove that these models all have the same free 
energy. Their thermodynamic behavior is determined by the spectrum of the transfer operator (Ruelle- 
Perron-Frobenius operator), which is defined using the maps (presentation functions) generating the Farey 
"tree" . The spectrum of this operator was completely determined by Prellberg. It follows that these models 
have a second-order phase transition with a specific heat divergence of the form C - [eln^ e]-^ . The spm 
chain models are also rigorously known to have a discontinuity in the magnetization at the phase transition. 
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I. INTRODUCTION 



In this work we consider several statistical models defined on the Farey fractions. One is the 
Farey fraction spin chain, a one-dimensional statistical model first proposed by two of the authors 
|||]. This work has spawned a number of further studies, by both physicists and number theorists 
|§, 1^ One can define the model as a periodic chain of sites with two possible spin states (A or B) 
at each site. The interactions are long-range, which allows a phase transition to exist in this one- 
dimensional system. The Farey spin chain is rigorously known to exhibit a single phase transition 
at temperature /3c = 2 Q. The phase transition itself is most unusual. The low temperature state 
is completely ordered ^] . In the limit of a long chain, for /? > /3c, the system is either all A 
or all B. Therefore the free energy is constant and the magnetization (defined via the difference 
in the number of spins in state A vs. those in state B) is completely saturated over this entire 
temperature range. Thus, even though the system has a phase transition at finite temperature, 
there are no thermal effects at all in the ordered state. 

At temperatures above the phase transition (for f3 < f3c), fluctuations occur, and the free energy 
decreases with /?. Here the system is paramagnetic. Since there is no symmetry-breaking field in 
the model, the magnetization vanishes. Thus the magnetization jumps from its saturated value 
in the low temperature phase, to zero in the high temperature phase [^]. This might suggest a 
first-order phase transition, but the behavior with temperature is different. In this work, we prove 
that as a function of temperature, the transition is second-order, and the same as that which occurs 
in the Knauf spin chain (see below) and the "Farey tree" multifractal model. The latter exhibits 
intermittency, and was studied by Feigenbaum, Procaccia, and Tel Q. 

The Farey spin chain is defined in an unusual way. It is given in terms of the energy of each 
possible configuration, rather than via a Hamiltonian. There is no known way to express the energy 
exactly in terms of the spin variables j^. Further, numerical results indicate that when one does, 
the Hamiltonian has all possible even interactions (and they are all ferromagnetic), so an explicit 
Hamiltonian representation, even if one could find it, would be exceedingly complicated. 

In previous work it was proven that the Farey spin chain free energy (per site, in the infinite 
chain limit) is the exactly same as the free energy of an earlier, related spin chain model due to 
Knauf 0. In the present paper, we extend this result in several ways. 

We begin by defining the spin chain and Farey tree models in Section II. In Section III we 
prove that the free energy for the Farey tree model is the same as the free energy of the Knauf 
model. This is established by use of bounds on the Knauf partition function. In Section IV, we 
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examine the Farey model, which is specified by the maps (presentation functions that generate 
the Farey tree. The free energy in this case is given by the logarithm of the largest eigenvalue 
A(/5) of the transfer operator. Some years ago, Knauf ||8| realized that the free energy of the Knauf 
model is also given by the logarithm of X{(3), (without noting the connection to the Farey tree 
model, however). Combining his result with our analysis rigorously shows the equality of all four 
free energies-for the Farey spin chain, Knauf model, Farey tree model and Farey model. In Section 
V, by using the results of we show that the phase transition is continuous (and of second order, 
i.e. the specific heat is divergent). It also follows that the phase transition in the Farey model 
occurs at the Hausdorff dimension of the Farey tree system, as expected. We conclude by briefly 
pointing out some connections with number theory and mentioning some implications of scaling 
theory for the spin chain models. 



II. DEFINITIONS 

We use the notation r^" := for the Farey fractions, where n is the order of the Farey fraction 
in level k. Level k = consists of the two fractions {fij}- Succeeding levels are generated by 
keeping all the fractions from level k in level k + 1, and including new fractions. The new fractions 
at level A; + 1 are defined via d^^^l := d^^^ + d^^'^^'^ and n^^^l := n["^ + n^^~^^\ so that 

^ = ill] 

k = '^ {t'^'t} 

k-? JO 1 1 2 u . 
1% — ^ \l'3>2'3'l/' 

Note that n = l,...,2'^ + l. When the Farey fractions are defined using matrices (spin states) A 
and B, the level k corresponds to the number of matrices and hence the length of the spin chain 

i- 

It follows that the fractions in a given level are always in increasing order. The Farey fractions 
differ from the Farey "tree" j^, where only the new fractions are kept at each succeeding level. 
The partition function for the Farey spin chain (FC) may be written as |||] 

Zf^(/3) := y \ , /3 G R. (1) 

Note from (||) that there are 2^ states at level k with energies -E^"^ = ln(d[."^ +n^"~''^''). The Farey 
fractions (and hence the energies) can also be defined using the spin variables A and B mentioned 
above W^, but this is not needed here. 



3 



For present purposes, it is convenient to use the partition function for the Knauf model 
which is rigorously known to have the same free energy as the Farey spin chain Q. The Knauf 
partition function may be defined via 

n=l )^ 

SO that a chain of length k has 2^^ states of energy = ln{d^f^^). The partition function can be 
written as sum of even and odd terms 

Zf(/3) = Z,^,(/3) + Z|f„(/3), (3) 

where 

2fc-l 2*~^ 

n=l l"fc ) n=l 1% / 

From the definition of the Farey fractions immediately follows 



and 



df-^)=da. (5) 



From (^) we have 



while from @ we obtain = Z^_-^{(3) so that 

Z,^,(/3) = Zf(/3)-Zf_i(/3). (6) 

The Farey tree model of Feigenbaum, Procaccia and Tel Q uses the "Farey tree" rather than 
the Farey fractions, which means retaining only the 2^~^ even fractions at level A; > 1 so we obtain 
the set 

{rf"V = l,---,2'="\fc>l}- 



The Farey tree partition function is defined by 



2fc-2 

} 

n=l 

The positive quantities (j'^f"''^ — r^^"" are the radii of the "balls" in this model. Note that we 
can also express this partition function using Farey tree denominators only. One finds 

9fc-2 / \ (5 

3 \ 



,(4n) ,(4n-2) 
n=l \"fc "fe 



III. EQUIVALENCE OF THE FAREY TREE AND KNAUF FREE ENERGIES 



In this section, we show the equivalence of the free energies of the Knauf and Farey tree models. 
We begin by finding bounds for the Farey tree partition function Z^{j3) in terms of the Knauf 
partition function. We are interested in the case /3 > 0, where there is a phase transition, but it 
will be easy to see that the free energies are equal for all /3 E M. 

The Farey fractions satisfy r^"^ = ^/('^i"'''^!" ^'')- This may be shown for instance using 

the matrix chain representation in Q. Thus 

(4n) _ (4n-2) _ (An) _ (An-l) (An-l) _ (An-2) 

r\j r\j r\j t\/ f\/ IX 

1 1 



1 

> 



4*"')' 



and similarly r^,"^"^ — r^"*" > 1/ ^4^^" ^"^ ) • Fi'o™ ® we also find 

(4n) {An-2) 2 



'(An-l)^^ 



Using (0) and (|), for /? > 0, gives 



fc-2 



ZkiP) > y — (10) 

p 2'=^^ / (An—2)\'^^ 

and also Zj^ > l2n=i ^/ [% j ■ Adding these two inequalities we find a lower bound for 



the Feigenbaum partition function 
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(") > I E j-^, = l^tm- (11) 

n=l [dl ' \ 



Using the inequality (pf) and the relation (0) gives the upper bound 



1 1 



Zf < 2^ y ^ = y ^ = l^Zti e(2/3). (12) 



Thus the Farey tree partition function at /? is bounded both above and below by the even part of 
the Knauf partition function at 2(3. 



^Z^^,{2P)<Z[{P)<2^Zti,em, (3>0 



(13) 
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Similarly, we can find, that 



2''zti,em < {(3) < ^<e(2/3), /3 < 0. (14) 



Finally, for /? = it is obvious that 



(/?) = Jzf (2/3). 



The free energy per site is defined by 



/(« JliM?). (15) 



(Recall that the level k corresponds to the length of the spin chain.) We now use (|T3|) to prove 
that 

fF{[3) = fKm- 

where fp refers to the free energy obtained from Z^. 
For P > 1 one has Q 

^K(^^. k^oo ({213 - 1) 

^^^^ ^ cm ' 

which implies that fx (2(3) = 0. Also, by (^), 

Z^J2(3)''^0, 



and using (13) gives 



Since Zf (/?) > 0, 



Z[{(3)'^0. 



k 

Note that for [3 = 1 one has Zf (1) < 1, since this partition function reduces to a simple sum of 
Farey tree fraction separations (ball lengths), which cannot exceed the length of the interval [0, 1]. 
Therefore the inequality still holds (and in fact, as shown below, /f(1) = 0). 
Now clearly 
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so by (O) we find 



and 



Thus we have 



1 1 



2(fc + l)2/3' 



-lnZ[{P) ^ 2/31n(fc + l) ^ ln2 
~ k k 



= /i^(2/5) = for /3 > 1. 



(16) 



The validity of /kC^) = is clear from the treatment in |^] and the remark at the end of this 
section. 



For P < 1 we can write 



SO 



k 



InZf ^^i^-i^' 



k 



k 



(17) 



It is shown in [|| (by arguments using the transfer operator, see below) that for < /? < 1 the free 
energies obtained from and Z^^ are the same, thus for k oo 



In (1 - ^ 



0. 



(18) 



(This also can be shown directly by considering the equation Z^{2f3) = 1 + '}2^j=i -^j,e^(2/?)) which 
follows from (^. For < /3 < 1 the series is bounded by a geometric series because of the inequality 
> 2^-^Z^_^^^.) For /3 < it is easy to check that Z^_^^^{2f3) / Z^^^{2f3) < 1/2. Thus (||) holds 
for all P < 1. 

Using (Il3|) (and, for /3 < 0, (111) and the line below) then establishes 



/f(/3) = fKm for P<1. 



(19) 



Note that, as mentioned, the Knauf partition function Z^{2[5) is finite as k ^ oo for /3 > 1 
Using (^) and ( |l3|) one sees immediately that the Farey tree partition function Z^{I3) vanishes in 
this limit for (3 > 1. At /? = 1, it follows immediately from the definition (^) and simple properties 
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of the Farey fractions that < Z[{1) < 1. For (3 < I, since /x(2/3) < [§ and using (|T|) and 
([l5| ) it follows that Z^{(3) is infinite. This establishes rigorously that the Hausdorff dimension of 
the set formed by the "balls" is Ph = 1, as expected. 

Finaly, consider ( [l3[) and the fact, mentioned above, that Z^{1) < 1. It follows that 

2^-1 

Zf,(2) = y —I— < 2, 

n=l j 

SO that this sum over the "new" Farey denominators is bounded by 2 at all levels. Since the "new" 
denominators at level k — 1 become "old" denominators at level k, one also sees that Zf (2) <2k+l. 

IV. TRANSFER OPERATOR APPROACH 

In this section we consider the transfer operator (Ruelle-Perron-Frobenius operator) of the Farey 
map. The previous section shows rigorously that the free energies of the Knauf and Farey fraction 
spin chain and Farey tree model are the same. Here we prove that they (as well as the free energy 
of the Farey tree model in a certain approximation specified below) are simply given by the largest 
eigenvalue of this operator. The next section considers the asymptotic behavior of this eigenvalue 



near the phase transition, known from the work of Prellberg |1C], which specifies the order of the 
phase transition. 

The Ruelle-Perron-Frobenius operator /C associated with a map / (piecewise monotonic trans- 
formation of closed interval /) is given by 

JCMx)= \ny)\-^ip{y), PeR, (20) 

where the sum is over each strictly monotonic and continuous piece of / satisfying the summation 
condition. See 0] for a more complete discussion. 
The Farey map is defined by ^ 

f Ux) = x/(l - x), < X < 1/2, 
[ /i(x) = (l-x)/x, l/2<x<l. 

The operator then consists of two corresponding terms /Cq and JCi which can be identified as 
"intermittent" and "chaotic" parts, respectively [|l^. We may write /C/3 = /C0 + /C1 where K,iip{x) = 
\Fl{x)\l^ <f{Fi{x)) and the "presentation function" Fi is the inverse map of fi (see (^) below). 
Thus 



K,,W = (l + .)-^[^(^)+^(^) 



/3 G R. (22) 
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Following the thermodynamic formalism approach [12| it was shown in |ri| that the largest 



eigenvalue of /C^ in ( ^21) (defined on the space of functions with bounded variation) is related to a 
free energy via f{(3) = —(3~^ In A(/3) for /? € M. We call this the free energy of the Farey model. 

In this section we consider )Cj3 acting on and show that the free energy obtained from its 
largest eigenvalue is the same as the free energy of the Knauf and Farey tree model (in its original 
version or using the approximation below) for < /? < 1. In the next section, we prove that the 
free energy of the Farey model in this /3 range is also the same. For (3 > 1, the free energy of any 
of these models is already known to be zero (see section III or ||lO[| ). 

The Knauf spin chain at level k — 1 may be described by a vector Yfe_i(2/3) G /^(Nq), the first 
component of which is the "even" Knauf partition function Z^^{2[3). The "transfer operator" of 
the Knauf spin chain then maps Yk-i{2(5) to the next level: 

yfc(2/5) =C(2/3)yfc„i(2/?), (23) 

where C(2/3) : /^(No) ^ ^{^0) and g 



-213 -i 



s=0 




(24) 



€ No), with the generalized binomial coefficients 



{Ii\zl{a-i))/b\, a G M, 6 G No, 



and 



if 6 < 0. Knauf § has further shown that for < /3 < 1, C{2(3) has the same 



largest eigenvalue A(/3) as : L^((0, 1)) L^((0, 1)). The argument involves expanding (22) 
about X = 1 with (/?(x) = X]m=o^™(-'^ ~ x)™. Doing this, one finds that the action of Kp on the 
quantities (note that am = (— 1)'"(^''™'^ (l)/m!) is given by C"^(2/3), where T denotes transpose. 

In addition, C'^{2(3) is independent of k, so the components of the vector Xk{2[3) (defined using 
(|2^ with C'^{2j3) replacing C(2/3)) are proportional to the Taylor series coefficients of an associated 



function (f)'^\x). This function therefore satisfies 



X) 



[l + x) 



-2/3 



X 



l + x 



+ 



l + x 



(25) 



It is shown in |] that C{2(3) (and hence C^(2/3)) is an operator of Perron- Frobenius type for < 
(3 < 1. Thus X{/3) is a simple eigenvalue (the same for C or C"^). The corresponding eigenvector is 
strictly positive and unique, and may be obtained (for C^) via V{2P) = limfc^oo ^ki'^P) /W^kC^P)]]- 
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In addition, it follows that for < /3 < 1 the eigenvalue A(/3) > 1 is an analytic function of /5, and 
its positive normalized eigenvector V{2f3) is analytic in p. Hence 

~ A(/3)V^^\ (26) 

where (f)^^\x) is the normalized eigenvector of ICp : L'^{{0, 1)) — > -^^^((0, 1)) corresponding to V{2P). 
Substituting this result in ( p5[ ) we get, for < /3 < 1, 

X \ ,IR\ ( 1 



A(/3)<^(^)(x) = + 



1 + X 



(27) 



which is equivalent to (22) when A(/3) is the maximal eigenvalue. Then 



together with dig), (|17|) and (|l^) give us the Knauf free energy as expected 



/;,(2/3) = -ilnA(/3), 0</3<l. (29) 

Note that for /? > 1, /x(2/3) = (see section III) and also that /(/3) = for /? > 1 follows 
from the spectrum of the operator ([0]> see also the next section). Thus the free energy of 
the Farey spin chain, Farey tree and Knauf models are given by the largest eigenvalue of the 
Ruelle-Perron-Frobenius operator for /? > 0. 

To further examine these connections we follow the treatment in |Q . We focus on ( p7D and make 
use of presentation functions. The Farey tree can be generated by two presentation functions 

1 - Fo = (30) 



1 + X 1 + X 

Every fraction at each level /c > 1 of the Farey tree can be reached by composition of k functions 
(e G {0, 1}) evaluated at x* = i. For example, at level /c = 3, Fq o Fi(i) = | = rg^^ So the 
diameter of every "ball" in the Farey tree model (see (^)) can be written as 



(4n) _ (4n-2) 
' k ' k 



o F,, o . . . o F,,_,(Fo(x*)) - F,, o F,, o . . . o F,,_,(Fi(x*))|. (31) 



Note that the sequence of presentation functions in the two Farey fractions in ( |3l|) is identical 
except for the F^j., i.e. only the presentation functions applied first to x* differ. As A; ^ oo, the 
diameter of the balls converges to zero (this follows easily from (|^). Therefore it is reasonable to 
suppose that for k sufficiently large each diameter can be approximated by the derivative of the 
composed function with respect to x*. Then, using the chain rule, (^) behaves asymptotically as 

^(4n) _ ^(4n-2) ^ j^,^ ^^^^ ^ ^ _ ^p,^ (F, O F,, O ...)... | . (32) 
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Thus we can write for the partition function 

~ • • • E ° ° • • oi^ E K^S^^^ ° ^^.+1 ° • • oi^ • • • (33) 

Notice that the sum over and ah lower indexed sums to its right depend only upon (F^^^^^ o 
Ftk+2 ° • • •)■ This motivates the definition 

V^f (x) := \Fl^{x)\P \PLM^.i^))f • • • ' (34) 
where (F^^. o F^^^^ o . . .) is denoted by x. One then finds 

(x) = Y \F'Ax)\Pi>f\{F,{x)). (35) 

e 

Note that since each presentation function i^^ is a ratio of polynomials, one can extend the definition 
of ^[.^''(x) to the whole interval [0, 1]. Substituting for F and F' we obtain ( [25| ) (with replacing 
Therefore choosing 'ip^^\x) > we find vi^^ ~^ V'''^'* as /c ^ oo, with the function ip^^'^ 
proportional to 0^^^ (the eigenfunction with the maximum eigenvalue A (/?)). This establishes that 
the approximation ( ^2] ) is exact in the limit k — > oo, as expected. 

Finally, it is interesting to note some connections with number theory. Specifically, for A = 1, 
(p7|) is known as the Lewis equation and has been studied (for complex /?) because of its connection 
to the Selberg (^-function and period polynomials (cusp forms of the modular group) [l^. An 
operator related to /C/j (^) also appears in this context and is called the Mayer operator [p!^]. 

V. ORDER OF THE PHASE TRANSITION AND DISCUSSION 

In the preceding, we have shown that the Farey spin chain the Knauf spin chain and 
(either version of) the Farey tree model all have the same free energy. Further, for < /? < 1 
their free energy is given by the largest eigenvalue of the Farey model transfer operator acting 
on (^)- Here we show that the transfer operator acting on the space of functions of bounded 
variation has the same leading eigenvalue in this /? range, which allows us to make use of the results 
of Prellberg. The corresponding equality of free energies for /3 > 1 (where the free energy vanishes) 
follows from known results, as remarked in the previous section. 

Prellberg has examined the spectrum of this operator acting on the space of functions with 
bounded variation |^ (details are in [|ll|). In order to make use of his results, we must show that 
the largest eigenvalue in this space is the same as that in L^((0, 1)). To prove this we examine 
the corresponding eigenvectors. Expanding ip{x), the eigenvector in the space, about x = 1 
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as above, one has ^{x) = J2m=o'^m{^ ~ x)'^. Thus ip{l) is finite, since the coefficients in this 
expansion are proportional to the components of the eigenvector of C'^{2(3) of largest eigenvalue 
(see section IV). Furthermore, the am are all positive, since the eigenvector of C'^ is positive. 
Therefore (p{x) is a (strictly) decreasing function on [0,1]. Finally, setting x = in ( p7[ ) shows 
that ip{0) is finite whenever A 7^ 1. Therefore, ip{x) is of bounded variation for < /3 < 1, and 
since both eigenvectors are unique (up to multiplicative constants) their eigenvalues must coincide 
in this range of /? values. 

The result of Prellberg of interest here is 

/?/(/?) = C|^^|^^ [1 + 0(1)], 0</3<l, 

where c > 0, and /?/(/?) = for (3 > 1. This form for the free energy is equivalent to that given in 
[^], as may be seen by use of the Lambert M/^- function. 

The non-analyticity at /? = 1 results in a phase transition of second order, since the second 
derivative of /(/?) diverges as (1 — /?)(ln(l — /?))^ as/3— >1~. This result agrees with where 
it is proven rigorously that the phase transition is at most second order. Note that the largest 
eigenvalue is discrete for /? < 1. For /? > 1, the discrete spectrum disappears and the largest 
eigenvalue becomes A = 1, which is the upper boundary of the continuous spectrum for all /?. 

Our result for the free energy also has some implications for the number of states of the spin 
chain models. The Knauf model partition function may be expressed as a Dirichlet series ^ 

00 

Z^{P) = J2Mn)n-^, (36) 

n=l 

where (j)kin) is non-zero when n is a Farey denominator at level k. This function converges from 
below to the Euler totient function <j){n) as /c — > 00. Since the energy of an allowed state is = In n, 
(pkin) gives the number of states of energy E at level k. The functions (pk and cj) are very irregular. 
Our result for the free energy then shows how the Dirichlet series in (^) diverges as /c ^ 00 for 
small (but positive!) (2-/3). (Recall that the phase transition in the spin chains appears at (3c = 2, 



since a factor of 2 appears in comparing with the Farey tree model, see (29).) For the Farey spin 
chain, an equation with the same form as ( |36|) may also be written, with the same leading divergent 
behavior. Here the limit of the function corresponding to (pki^) is not known, though some related 
information is available [^. 

One can also consider the implications of scaling theory for the two spin chain models. It is 
known that the magnetization (defined via the difference in the number of spins in state A vs. 
those in state B) is one for temperatures below the transition and zero above it |, ||. Thus the 
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magnetization jumps from its fully saturated value to zero at the transition. This would lead one 
to suspect a first-order transition, but as we have seen, the behavior with temperature is second- 
order. However, both these results seem to be consistent with scaling theory, with renormalization 
group eigenvalues yx = d and yn = d, where d is the dimensionality, and using (2 — /5)/ln(2 — (3) 
as the temperature scaling variable. We plan to report more fully on this elsewhere. 
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